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ABSTRACT 

We present Keck/DEIMOS spectroscopy of resolved stars in the M31 satellites And XXVIII & 
And XXIX. We show that these are likely self-bound galaxies based on 18 and 24 members in 
And XXVIII & And XXIX, respectively. And XXVIII has a systemic velocity of -331.1 ± 1.8 km s'^ 
and velocity dispersion of 4.9 ± 1.6 km s~^, implying a mass-to-light ratio (within ri/2) of ~ 44 ± 41. 
And XXIX has a systemic velocity of — 194.4± 1.5 km s^^ and velocity dispersion of 5.7± 1.2 km s~^, 
implying a mass-to-light ratio (within ri/2) of ^ 124±72. The internal kinematics and implied masses 
of And XXVIII & And XXIX are similar to dwarf spheroidals (dSphs) of comparable luminosities, 
implying that these objects are dark matter-dominated dwarf galaxies. Despite the large projected 
distances from their host (380 and 188 kpc), the kinematics of these dSph suggest that they are bound 
M31 satellites. 

Subject headings: Local Group — galaxies: dwarf — galaxies: individual (And XXVIII, And XXIX, 
M31) — galaxies: kinematics and dynamics 



L INTRODUCTION 

Local Group (LG) dwarf galaxies are crucial el- 
ements of near-field cosmology. They serve as 
valuable probes of low luminosity galaxy formation, 
and provide an important windo w into the nature 
of satellite- host interactions (e.g.. Bullock ct al. 2000'; 
Strig ar i ct al- .2007: .Kravtsoy. ,2010 : Ka zantzidis ct al. i 
2011t iBovlan-Kolchin et all 120111 : lAjod crhalde n et al.l 
20121 ) . Their utility stems primarily from their proximity 



and low surface densities, allowing spectroscopic observa- 
tions of individual stars. This enables far more in-depth 
studies than are possible for more distant targets. 

The numbers of known M31 and Milky Way (MW) 
dSph satellites has increased greatly in recent years. This 
is predominantly due to homogenous wide-field surveys, 
particularl y the Sloan Digital Sky Survey (SDSS) fo r 
the MW dWillman et al.l ^OOSt IBelokurov et all [2007I) . 
and the Pan- And romeda A rcheological Surve y (PAndAS, 
llbata et all 120071 : iMcConnachie etaH l2009( ) . The envi- 
rons of M31 have been particularly fertile ground for dis- 
covering dSphs, yielding 20 of its 28 known dSphs in the 
last six years. 

Two of the most recently discovered M 31 satell i te can - 
didates are And X XVIII, found by Slat er et al.l ()2011[ ). 
& And XXIX, bv iBell et all (f201ll ). Projected in the 
plane of the sky at the distance of M31, And XXVIII & 
And XXIX lie at 380 and 188 kpc, respectively. The 3D 
distances from M31 are 367 and 188 kpc, using the tip 
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of the red giant bra nch (TRGB) line-of-si ght distances of 
iSlater et al.l (|20T1 and IBell et all poTll ). The expected 
virial radius for M31's dark matter hal o is ~ 300 kpc 
(|KlvT)in et al.l[200l Imtkins et aLllMoh . so And XXIX 
lies in the outskirts of M31's halo and And XXVIII just 
outside. This makes And XXVIII the second-most dis- 
tant M31 dSpl0 satellite, with only And XVIII more 
distant (see !j3.3|) . 

And XXVIII, in particular, is potentially an impor- 
tant data point for galaxy formation. It may be an 
analog of Leo T, a satellite in the outskirts of the MW 
halo and the only low- luminosity MW satellite known 
to co ntain HI gas lllrwin et"aLl l2007t iRvan- Weber et all 
120081 ). This motivates more detailed studies of these 
M31 satellites, as Leo T has been valuable for under- 
standing star formation and gas stripping in dSphs (e.g., 
Ricottill2009HGrcevich fc Putman 2009; B ovill fc Ricottil 



20111: INichols fc Bland-Hawthorn.,2011.) . If they have no 



recent star formation. And XXVIII & And XXIX may 
be more akin to Tucana and Cetus. These are unusual, 
non-starforming dS phs of the LG that do not seem to be 
satellitesjLavcrv fc Mighelll 119921: IWhiting et all 119991: 
IMonelh etal.. .20101) . but may be "backsplash" galaxies 
that were once much closer to either M31 or the MW 
(|Knebe et al.ll20nl: lO^n et al.ll2013l) . 

Here we present the first spectroscopic observations 
of And XXVIII & And XXIX, obtaining kinematics of 
resolved stars using the DEIMOS spectrograph on the 
Keck II Telescope. In 321 we describe our observations 
and analysis procedures, in Jj3]we present our kinematical 
measurements for And XXVIII & And XXIX, and in 311 
we summarize and conclude. Throughout this paper we 
adopt a distance modulus to M31 o f mm.si = 24.47, cor- 
responding to a distance of 783 kpc (|McConnachie et aLl 
[inni), and distances of 650^80° (a* = 24.06t;^-^) and 
730 ± 75 kpc = 24.31 ± 0.22) for And XXVIII & 

^ The starforming dirrs IC 1613 and Peg are more distant and 
possibly associated with M31I IJVlcConnachic 2012). 
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And XXIX, respectively ([Slater et all 1201 H iBell et alj 

[gonl). 

2. OBSERVATIONS AND ANALYSIS 

2.1. Target Selection and Reduction 

We selected stars for spectroscopic follo w-up from the 
SPSS Data Release 8 photometric catalogs ()Aihara et al.l 
[20T1 . We selected stars near And XXVIII & And XXIX, 
eliminating objects where the SDSS model and psf mag- 
nitudes were discrepant by more than 0.25 mags. We 
then selected candidate dSph stars near the red giant 
branch (RGB) of a fiducial isochrone in the g~i, r color- 
magnitude diagram (CMP). W e used 12 Gyr Partmouth 
isochrones (jPotter et al.l I2008D offset by ISchlegel et al.l 
(|1998f ) extinctions with [Fe/H] and distance m odulus 
matched to Arid XX VIII & And XXIX (,Slater et al.l 
120111: iBeh et al.l 120 111) . We then populated the slitmask 
by selecting stars in the following sequence, prioritizing 
brighter stars for each group: those within 0.3 mags of 
the isochrone, those within 0.6 mags, and all remaining 
stars. The resulting selections covered a wide enough 
area in the CMP that our results do not dep end strongly 
on the ass umed dista i ices a nd [Fe/H] from iSlater et al.l 
((20T1 and lBell et aT\ ([20lTh . 

2.2. Observations and Reduction 

Spectroscopic observations were obtained on the nights 
of April 22-23 and September 16-17, 2012 using 
the PEIMOS spe ctrograph on the Keck II telescope 
(|Faber et al.ll2003[) . We used the 1200 lines mm"! grat- 
ing covering a wavelength region of 6400 - 9100 A. This 
provided a FWHM resolution of w 1.3 A, equivalent to 
50 km s~^ at the center of our wavelength range. We 
observed 3 slitmasks for And XXVIII, and 2 slitmasks 
for And XXIX, with an average exposure time of 3200 
sec per mask. 

Our spectroscopi c reductions cl o sely f ollow the pro- 
cedure outlined in iTollerud et al.l (j2012l) . We reduce 
the spectra using the spec2d pipeline developed for 
the PEEP2 sur vey, which produces Id spectra from 
the raw image s (iPavis et all [200l iCooper et aL|[20T2l: 
INewman et al.l I2012I ). In Figure [TJ we show smoothed 
co-adds of the resulting spectra for members stars of 
And XXVIII & And XXIX (see g23]and|3]for member- 
ship determination). The average signal-to- noise ratio 
for these spectra was 4.8 per pixel. 

We cross-correlate these Id spectra with a series of high 
signal-to-noise ratio (SNR) templates of known radial ve- 
locity to determine a line-of-sight velocity. We determine 
errors on these velocities by re-simulating each spectrum 
1000 times with noise added to simulate the per-pixel 
variance. We measure a cross-correlation velocity for 
each re-simulation, and use the mean and standard devi- 
ation of the distribution of velocities as the line-of-sight 
veloc ity and uncertainty for each star (jSimon fc Gehal 
I2007D . We also add in quadrature a systematic floor to 
the uncertainty, determined by repeat measurements of 
star s from Sim on k Geha (2007), Kalirai ct al. (2010), 
and ITollerud et al.l (|2012D . These each found consistent 
values of of 2.2 km s~^ for the systematic floor over a 
temporal baseline longer than th e time from these obser- 
vations to ITollerud et all (|2012D . Hence, we apply this 
same floor to the measurements in this work. 
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Figure 1. Continuum- normalized spectra of And XXVIII &: 
And XXIX member stars. These spectra are inverse vari- 
ance weighted hehocentric frame co-adds of 18 and 24 stars for 
And XXVIII & And XXIX, respectively, and have been smoothed 
with a 3-pixel gaussian filter for clarity. The lower (blue) line is for 
And XXVIII, while the upper (green) line is And XXIX. The top 
panel is a zoom near Ha, while the lowe r pan el is a window that 
includes the calcium triplet (CaT). See i|2.3l and l3l for details of 
membership determination. The vertical black lines marking Ha 
and the CaT are at wavelengths where the features would be for 
a velocity halfway between the Dsys of And XXVIII & And XXIX. 
The offsets of the features in the spectra from these wavelengths 
are due to the differing Vaya of the dSphs. 



2.3. Membership 

A clean sample of member stars is crucial for determin- 
ing the internal kinematics of And XXVIII & And XXIX. 
The two main sources of contamination for these obser- 
vations are foreground MW stars and M31 halo stars. 
As demonstrated below, the MW foreground star veloc- 
ity distribution is disjoint from the dSphs in the direc- 
tion of And XXVIII & And XXIX, and hence the MW 
contaminants can mostly be filtered with velocity cuts. 
And XXVIII & And XXIX are far from M31 in projec- 
tion (^ 380 and ^ 190 kpc projected, respectively), and 
thus the M31 halo dominates over the bulge or disk com- 
ponent (jCourteau et all 120111: iGilbert et all I2012D . We 
demonstrate in W2A\ that even the M31 halo density is 
low enough that its contamination is unlikely to affect 
our kinematic results for And XXVIII & And XXIX. 

For both And XXVIII & And XXIX, our membership 
determination begins by selecting stars in the r — i, r 
CMP. We select stars in a CM P box that encompa sses 
stars near PARSEC isochrones (jBressan "eFaI1l20ll) for 
a 12 Gyr populatio n with [ Fe/H] values that match the 
estimates from , Slater et al.l (f2011, ) and .Bell et al.. (,2011i) . 
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Table 1 

Photometric and Kinematic Properties of And 28 & And 29 



Row Quantity Units And XXVIII And XXIX 



(1) 
(2) 
(3) 

(4) 
(5) 
(6) 
(7) 

(8) 



a (J2000) 
<5 (J2000) 
^LOS 
Mv 

RcH 

RcS 

[Fe/H] 
'1/2 



n : m : s 
o . / . // 

kpc 

mag 

arcmin 

pc 

dex 

pc 



22''32™41!2 
31° 12' 58.2" 
65011^° 

1.11 ± 0.21 
210+60 

2.0 



23'^58™55?6 
30°45'20.0" 
730 ± 75 
-8.3 ± 0.4 
1.7 ± 0.2 
360 ± 60 
~ -1.8 
480 ± 80 





Keck/DEIMOS Results 




(9) 


-^member 18 


24 


(10) 


Vsys kms-1 -331.1 ±1.8 


-194.4 ± 1.5 


(11) 


o-LOS kms-1 4.9 ±1.6 


5.7 ± 1.2 


(12) 


M/Lv(<n/2) Mq/Lq 44 ±41 


124 ± 72 


(13) 


Afi/2 Mq X 10** 4.7 ±3.2 


11.1 ±4.9 



Note. — (1-7) Right ascension, dcclinatic 
tnde, effective radii, and metallicity are talte 
iiiOlii ) for And 28 & And 29, respectively. 

These boxes are shown in the upper-left panels of Figures 
[2] and [21 We place the upper edge of these boxes above 
the tip of the isochrone because this accepts TRGB stars 
even if the And XXVIII & And XXIX distance moduli 
are ~ la off. To filter out the MW foreground stars that 
lie near the dSph members in the CMD, we impose a ve- 
locity window centered on the "cold spikes" apparent in 
the lower- left panels of Figures [Hand [21 The stars within 
these peaks and inside the CMD box are spatially con- 
centrated near the centers of And XXVIII & And XXIX, 
as expected for a self-bound galaxy. Additionally, the 
absence of stars near these velocities above the TRGB 
of our chosen isochrones show even 2a errors in the dis- 
tances to And XXVIII & And XXIX would have no effect 
on our membership determination. 

The velocity distributions of these candidate members 
are consistent w it h a G aussi an in the sense that b o th the 
IShaoiro & Wiig (fl965l ) and lAnderson fc Darlind (fl952l ) 
tests cannot reject the null hypothesis of Gaussianity at 
the p — 0.05 level. While these Gaussian peaks imply 
that most of the selected stars are members, we cannot 
discount the possibility of a small number of contami- 
nants that by chance are near the dSph locus in both 
the CMD and velocity. To address this, we estimate the 
surface density of MW foreground stars expected in our 
spectrosc opic sampl e using the Besangon model of the 
MW (jRobin et all 12003). We select model stars in the 
direction of each dSph which fall within the range of 
our observations in the r — i, r CMD. We then normal- 
ize the resulting model population such that the overall 
number of stars in the model match our observations for 
"yiiciio > —100 km s^^, far from the dSph velocity. Fi- 
nally, we determine how many stars in the model would 
fall within the velocity window for the dSph members. 
For both And XXVIII & And XXIX this number is small, 
~ 0.7 and ~ 0.8, respectively, in the entire field of our 
observations. 

2.4. Kinematical Parameters 

With a member sample selected as described above, 
we proceed to determine the internal kinematics of 
And XXVIII & And XXIX. In the sections below, we 



c-of-siglit distance, absolnte magni- 
a lSiater et al.| j^Olij ) and lijeil et al.| 



model the velocity distribution of the member stars as 
a Gaussian with mean Vsys (systemic velocity), and dis- 
persion iJlos • We determine the kinematical parameters 
of this model with a maximum likelihood fit, weighting 
all members equally, and use the inverse of the Hessian 
matr ix to estimate the uncertainties on thes e parame- 
ters f Walker et all [20071 : I^Hlerud et al.l[20Tl) . We ex- 
amine the resulting likelihood maps for Vsys and ctlos, 
and note that they are very close to Gaussian within la 
of the peak. This validates the use of the Hessian matrix 
for estimating uncertainties, as it implicitly assumes the 
likelihood function near the maximum likelihood peak is 
approximately normal. 

This approach assumes each star is only a tracer of 
the internal kinematics of the galaxy, and hence does 
not acco unt for the i mpac t of binary stars. As de- 
tailed in [Minor et all (|2010( ) and iMcConnachie fc Cotel 
(|201Clf l. unresolved binaries can inflate the internal ve- 
locity dispersions measured for dSphs. Repeat veloc- 
ity measurements are required to provide any meaning- 
ful correction for this effect, and this is not available 
for our data set. However, the uncertainties we report 
for And XXVIII & And XXIX below are larger than 
the corrections that result fro m multi-epo ch observations 
of MW dSphs (Mi nor et all 2010; Marti nez etaD 120111 : 
[Minor! '2013). Hence, unless M31 dSphs have substan- 
tially different binary populations than MW dSphs, is it 
unlikely unresolved binaries would have a major effect on 
the results we present here. 

We also estimate the M31 halo surface density near 
And XXVIII & And XXIX. We start by determining 
stellar luminosity functions from PARSEC isochrones 
for [Fc/H] = —0.5 and —1.5, a range that approxi- 
mately samples the metallicity range of the M31 halo 
([Gilbert et al.l[20l"^ . We use these luminosity functions 
to determine how many stars bright enough to meet our 
spectroscopic survey limits (r < 22.5) are present for a 
given surface brightness. For each dSph we then deter- 
mine the M31 h alo surface br i ghtnes s using the best- fit 
power law from [Gilbert et aTl ([2012D . Combining these 
surface brightness estimates with the luminosity func- 
tion thus provides an estimate of the number of M31 
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Figure 2. Upper-left: Foreground extinction-corrected r, r — i CMD of the SDSS photometry of stars in the And XXVIII field. The 
points (black) are all star-like objects in the SDSS catalog in the spatial and CMD regions shown. Stars with spectroscopic data are co lored 
by t^helioi with solid circles representing stars that are classified as members, outlined triangles are possible members (see text in i|3.1l l. 
and asterisks are non-members. The (black) dashed box indicates the CMD selection window for membership, and the black error bars 
at r — i = 1.5 are photometric errors in the SDSS averaged in magnitude bins. The solid (red) line is a PARSEC isochrone of 12 Gyr 
age and [Fe/H]= —1.5 (Bressan et al. 2012). Upper-right: Spatial Distribution of stars in the And XXVIII field. Symbols here are the 
same as those in the upper-left panel. The solid (red) line indicates the half-light ellipse. Lower-left: Velocity histogram for stars in 
the And XXVIII field. The blue histogram is member stars, while red is f or non-members. T he (black) dashed vertical lines indicate the 
velocity window used for membership. The gray line is the histogram of the lRobin et al.l II2003I ) model of MW foregrounds stars, normalized 
as described in i|2.3l Lower-right: Radial velocity vs. distance from dSph in the And XXVIII field. Filled circles (blue) are members, 
outlined triangles (yellow) are the two stars of uncertain status, and star-shaped symbols (red) are non-members. Solid vertical lines are 
per-star velocity uncertainties. The (black) dashed line is the velocity window used for determining membership, while the (red) dotted 
line indicates the half-light radius of the dSph. 

halo stars expected in our fields. We estimate ~ 0.3 
M31 halo stars in our entire And XXVIII field, and ~ 2 
M31 halo stars in our And XXIX observations. Further- 
more, the M31 halo velocity distribution, while overlap- 
ping with the dSphs, is much hotter (cr ~ 100 km s~^, 
IChapman et al.l I2OO60 . further reducing the likelihood 
of an M31 halo star lying within our velocity window. 
While the predicted M31 halo and MW foreground con- 
tamination is small, we nonetheless simulate its effects 
.41 showing that it should have no impact on our 



inferred kinematical results. 



3. KINEMATIC RESULTS 
3.1. And XXVIII 

We present our observations of stars in the 
And XXVIII field in Figure [5] This object is far from 
M31 (~ 380 kpc projected), but is relatively close to the 
Galactic plane {b ~ — 23deg). There is a clear velocity 
peak at Ugys ~ —300 km s~^ (lower-right panel of Fig- 
ure [2]), well away from the velocity peak for MW stars 
(^hoiio ~ 0). Applying our membership and kinematical 
analysis described in ^ 32. 31 and 12.41 (using a velocity win- 
dow of -350 to -300 km s^i) yields Vsys = -328.0 ± 2.3 
km s~^ and ctlos — 8.1 ± 1.8 km s~^. Member stars are 
clustered near the photometric center (Upper-right panel 
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Figure 3. Same as Figure [2] but for And XXIX. 

of Figure!!]), and generally lie close to a fidueial isochrone 
for a dSpli-like stellar population (Upper-right panel of 
Figure[2]). While there are a few stars with velocities near 
the cold peak that appear to be outliers in the CMD, this 
is likely due to the relatively high photometric errors in 
this SDSS field. 

The right panels of Figure [5] reveal two stars catego- 
rized as members that are at large distances from the 
photometric center of And XXVIII. Based on the analysis 
desc ribed in ^2. 31 (assum ing the surface brightness profile 
from'Slate r et al.ll201lD . there should be no And XXVIII 
RGB stars at projected distances farther than these two 
stars (~ 1100 pc). These two stars are also velocity out- 
liers. Thus, either And XXVIII's surface brightness pro- 
file is intrinsically very different from other dSphs, the 
stars were tidally stripped, or these two stars are not as- 
sociated with And XXVIII. Both stars' spectra have an 
absorption equivalent width of ^ 0.7 A for Na I A8190, 
a surface-gravity sensitive feature that is only detected 
in dwarf stars. This equivalent width marginally sug- 
gests they are dwarf stars, and hence unassociated fore- 
grouncfl. While the foreground star estimate from H'l.'Sl 

^ For one of these stars tlie measurement may be contaminated 
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suggests that two MW foregrounds are unlikely, Poisson 
statistics imply a probability of ^ 5%, so it is possible 
that both stars are contaminants. 

Hence, we also determine the kinematical parameters if 
these two stars are excluded, yielding Vsys — —331.1 ±1.8 
km s""'^ and clos — 4.9 ± 1.6 km s~^. ctlos in this case 
is significantly smaller than if these stars are included, 
but is stable within Icr to the exclusion of any other 
two stars. In the hypothesis that these stars are not 
members, this velocity dispersion provides an estimate 
of And XXVIII's mass (see S|33] for details). We com- 
bine this mass with the 3D distance from M31 (367 kpc) 
and a mass for M31 of 2 X (intentionally in the 

high range of M31 mas ses) to estimate the Jacob i tidal 
radius of And XXVIII f Binnev fc Tremaind 120081 ). 

We find a tidal radius for And XXVIII of ~ 3.4 kpc, an 
order of magnitude larger than the current stellar extent 
{Re ~ 0.21 kpc). This strongly supports the hypothe- 
sis that the two outlier stars were not recently stripped. 
Further, the Vgys measured here for And XXVIII is very 
close to the M31 Wgys- This may imply that it is near 
apocenter in its orbit, and thus is moving slowly (e.g., 

by a sky line. 
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iBovlan-Kolchin et alj I2012D . Hence, stars stripped at 
pericenter should be far from And XXVIII. Taken to- 
gether, this suggests And XXVIII has not recently been 
tidally stripped, and the two anomalous stars are con- 
taminants. For our fiducial kinematic measurements of 
And XXVIII (listed in Table [1]), we exclude these two 
stars as likely contaminants. 

In either described below in i )3.5[ this dis- 

persion implies a mass-to-light ratio larger than any 
plausible stellar population. Hence, we confirm that 
And XXVIII is most likely a self-bound galaxy with a 
massive halo. Its status as a satellite of MSI is examined 
in detail in ^'S.'3^ 

3.2. And XXIX 

We present our spectroscopic measurements of stars in 
the And XXIX field along with the corresponding SDSS 
photometry in Figure [31 As in ^13. 11 a clear cold peak 
is present, but for And XXIX, the peak is closer to the 

MW peak at WhoUo 200 km s^^. While And XXIX is 

closer to M31 than And XXVIII 190 kpc projected), 
it is more offset from M31's Wgys, and is farther from the 

Galactic plane than And XXVIII {b 31deg). Thus 

both the MW and M31 contamination for the And XXIX 
field is less than And XXVIII. Indeed, applying our mem- 
bership methodology to the cold peak (using a veloc- 
ity window of —215 to —175 km s~^) provides the his- 
togram shown in the lower-right panel of Figure [3l This 
is consistent with a Gaussian distribution and has no ob- 
vious velocity outliers at large distances. We m easure 
its kinematical parameters as described in ^2A\ deter- 
mine Usys = — 194.4 ± 1.5 km s~^ and ctlos — 5.7 ± 1.2 
km s~^, and tabulate those results in Table [1] As for 
And XXVIII, these kinematical parameters are similar 
to those of other M31 and MW dSphs. We conclude 
that And XXVIII is also a self-bound dSph satellite of 
M31. 

3.3. M31 Satellites or Local Group Field? 

Both And XXVIII & And XXIX are at relatively large 
projected distances from M31 of 380 and 190 kpc, re- 
spectively. By contrast, the majority of M31's dSphs lie 
within 150 kpc projected. This large difference is primar- 
ily a selection effect caused by the angular extent of the 
PAndAS survey. Nevertheless, the large projected sepa- 
ration begs the question of whether or not And XXVIII & 
And XXIX are truly satellites of M31 rather than "free 
fioating" galaxies of the LG (e.g., Tucana and Cetus). 
While the typical usage of the word "satellite" for these 
galaxies is sometimes ambiguous, here we define it to 
mean an object that is instantaneously bound to its host 
at the current epoch. 

Figure |4] plots Vgys, the line-of-sight relative to M31, 
against the 3D distance from M31 for a selection of 
dSphs. Also shown (as green curves) are the one- 
dimensional esc ape velocity ( i. e., fps r. /VS) for the M31 
halo model of iKlvpin et ahl (|2002[ ) (Mvir = 1.6 x 
IO^^Mq). a satellite lying above this curve is thus un- 
bound if its two tangential velocity components are equal 
to its observed line-of-sight Vgys- It is clear from inspec- 
tion of this plot that both And XXVIII (blue circle) & 
And XXIX (red square) lie well within these curves, in- 
dicating that their tangential velocity components must 




£?3D,M31 [kpc] 



Figure 4. Distance vs. line-of-sight Vb-vb for M31-associated 
dwarf galaxies. The (blue) circle and (red) square are from this 
work, And XXVIII & And XXIX, respectiv ely. Up-point i ng tr i- 
angles (gray) are dSphs fr om the samples of ICoUins et al.l 11201 II ), 
iToUeru d et al. (2012) and IHo et all ijoig V Diamo nds and right- 
pointi ng triangles (gray) are dirrs and dEs from IMcConnachid 
1120121 ). The distance axis is 3D distance from M31, computed by 
combining the projected separation of the dSphs from M31 with 
line-of-sight TRGB distances. T he (black) vertical dotted line is 
the virial radius of M31 from the'Kly pin e t al. (2002) model. The 
(green) curves are the one-dimensional escape veloci ties f r om a n 
NFW dark matter halo matching the model o fKlvpin et al.l II2002I '). 
To make them one-dimensional, we divide the escape velocities by 
\/3, so they represent the Vbjb above which a satellite would be 
unbound if its (unmeasured) tangential velocity components were 
equal to jJsys . And XXVIII & And XXIX lie well within the bound 
region. 

exceed the line-of-sight velocity by a large margin to not 
be bound to M31. This holds even for And XXVIII, 
despite the fact that it lies beyond M31's virial radius 
(vertical dashed line in Figure 31 

Also of note in Figure 3] is the satellite status of the 
most distant dSph, And XVIII. Its Wsys is only ^ 30 
km s^^ from M31, despite being nearly as far behind M31 
as M31 is from the MW. While this strongly suggests it is 
formally bound to M31, the free-fall time for And XVIII 
to reach M31 is greater tha n the time in which M31 i s 
likely to merge with the MW ()van der Marel et al.ll2012D . 
Hence, systems as distant from M31 (or the MW) as And 
XVIII might be better considered satellites of the future 
merged M31/MW galaxy. 

While the Wgys measurements for these galaxies are 
suggestive, they cannot definitively determine if a satel- 
lite is bound due to the unmeasured tangential veloc- 
ity. However, the evidence that dark matter halos are 
present in both And XXV III fc And XXIX ( MM pro- 
vides strong constraints. IBovlan-Kolchin et a l. (20121) 
shows that in ACDM simulations of either MW or M31- 
like halos, halos at distances from the host like those of 
And XXVIII or And XXIX are nearly always (> 99.9 %) 
bound. In a ACDM context, it is thus extremely likely 
that And XXVIII & And XXIX are M31 satellites. 
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3.4. Contamination Simulations 



The method described in t j2.4l for estimating kinemat- 
ical parameters assumes a pure sample of dSph member 
stars and does not formally account for the possibility 
of contamination due to MW foregrounds and M31 halo 
stars (discussed in i j2.3[) . To determine if such contami- 
nation might affect our kinematical results, we simulate 
the effect of contamination on our parameter estimates 
by creating mock velocity datasets designed to mimic 
And XXVIII & And XXIX. 

For each dSph we create 10000 mock datasets com- 
posed of a Gaussian distribution with kinematic param- 
eters corresponding to our results for the dSph. We add 
to this two additional uniform distributions with nor- 
malization set to match the number of expected stars in 
the MW foreground and M31 halo as estimated in §2.31 
We account for fractional numbers of contaminants by 
only including an additional contaminant star for a cor- 
responding fraction of the mock datasets. We then mea- 
sure the kinematical parameters of these mock velocity 
distributions, and compare them to the input parameters 
used to generate the mock datasets. 

We generate mock datasets based on both the 
And XXVIII & And XXIX fields following the prescrip- 
tion described above. For both sets of mocks, the esti- 
mated (Tlos and Wgys were very close to the true value 
used to initialize the mock datasets. More specifically, 
the variance in the difference between the true and esti- 
mated values due to the contamination is ^ 5x lower than 
the variance in each measurement due to small number 
statistics. Thus, contamination at the level we estimate 
here cannot have a statistically significant impact on our 
results. 

3.5. Mass Estimates 

With reliable kinematical parameters in hand, we are 
now in a position to estimate the mass of these galaxies 
from their internal velocities. We search for rotation by 
dividing the galaxies in half and search for mean veloci- 
ties of opposite sign on each side, over a range of position 
angles. We find that the maximal rotation signal for ei- 
ther galaxy is negligible relative to the dispersion, < 1 
km s~^. Hence we model the galaxies as purely pressure- 
supported systems. We estimate the mass of these 
galaxies from their velocity dispersions {M1/2) following 
Equation 2 of iWolf et al.l ()2010D . and tabulate these for 
And XXVIII & And XXIX in Table [H The mass ob- 
tained with this formula is not strongly degenerate with 
the anisotropy, but in this interpretation, it is only valid 
as a mass measurement within the deprojected (3D) half- 
light radiufQ. Applying this to our kinematical parame- 
ter estimates yields Mi/2(< ri/2) = 4.7 ± 3.2 x IO^Mq 
for And XXVIII, and Mi/2(< r^/^) = 1.1±0.49 x IO^Mq 
for And XXIX. 

A simple mass-to-light ratio estimate may be obtained 
by dividing this M1/2 by the luminosity within ri/2 ~ 
by definition this is L/2. For And XXVIII, we obtain 
M/Lv{< ri/2) = 44±41 Mq/Lq, and for And XXIX, we 
find M/Lv{< ri/2) = 123 ± 72 Mq/Lq. Like other M31 

''' For surface brightness profiles appropriate for dSphs, however, 
the deprojected (3D) half-Hght radius is well estimated as 4i?(,ft/3 
where i?(,ft is the projected half-light radius. 
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Figure 5. Log of Mass-to-light ratio vs. Luminosity for M31 
dSphs. And XXVIII & And XXIX are represented as a (blue) cir- 
cle an d (red) square, while other M31 dSphs from Fl^Uerud et al.l 
II2OI2I) and[Ho et al. (2012) arc (gray) triangles. Masses are deter- 
mined using the Wolf et al. (2010) mass estimator (using half-light 
radii from 'Brasscur et al. 2011), which is valid when interpreted 
as the mass within the deprojected (3D) half-light radius. Both 
And XXVIII & And XXIX lie in the range of other M31 dSphs and 
have mass-to-light ratios greater than plausible from their stellar 
populations alone, implying they are dark matter-dominated. 

dSphs of similar luminosity (see Figure [5]), And XXVIII 
& And XXIX have mass-to-light ratios above those that 
are possible for purely stellar systems. This reveals 
the presence of a dark matter halo, even within the 
half-light rad ius, where baryons d o minate for brighter 
galaxies (e.g.. 



iToUerud et al 



Strigari et al1l2008bllal : IWalker et al.l[200l 
■|2011D. For And XXVIII, the uncertainty 



admits a stellar-like M/L at the la level, primarily due 
to uncertainty in the galaxy's luminosity, but the best 
estimate is well above any plaus i ble st ellar value. Al- 
ternatively, [McGauiF&Miigro^ ()2013[ ). under a Mod- 
ified Newtonian Dynamics (MOND) hypothesis, have re- 
cently predicted velocity dispersions for And XXVIII & 
And XXIX of 4.3t5];? ^■'^-0.7 km s'^ (error bars are 
from assuming mass-to- light ratios from 1 to 4 Mq/Lq). 
These predictions are marginally consistent with our 
measurements, but whether this is numerical coincidence 
given the uncertainties or evidence for MOND is be- 
yond the scope of this paper (Kap linghat & Turner 200j; 
IMcGaughll20ni : [Foreman fc Scottll2012t lGnedinil2012f) . 

In Figure we show the mass-to-light ratios of 
And XXVIII & And XXIX in the context of other M31 
dSphs. This plot demonstrates that their internal dy- 
namics are fairly typical of M31's dSphs, despite their 
distance from their host. This implies that the processes 
shaping these dSphs' kinematic structure either operate 
quickly upon infall, are in place before the satellites inter- 
acted with M31, or both And XXVIII & And XXIX have 
been in the M31 system for some time and have already 
been influenced by M31's environment fe.g.. lWetzel et aLl 
l2012f) . These trends are also similar to MW dSphs, sug- 
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Resting these process es are not unique to the M31 system 
(|Tollerud et aHIMl . 

4. CONCLUSIONS 

The spectroscopic observations of the And XXVIII & 
And XXIX fields described here support the foUowing 
conclusions: 

1. We have spectroscopically confirmed that 
And XXVIII & And XXIX are likely self- 
bound dwarf galaxies, with velocity dispersions of 
4.9 ± 1.6 and 5.7 ± 1.2 km s"\ respectively. The 
implied large mass-to-light ratios are consistent 
with dark matter-dominated dynamics. 

2. While they are in the outskirts of the M31 sys- 
tem, the systemic velocities of And XXVIII & 
And XXIX imply they are both bound to M31. 

3. The internal kinematics and implied masses of 
And XXVIII & And XXIX are like those of other 
dSphs of similar luminosities, despite the large dis- 
tances from their host. 

While these results clearly demonstrate that 
And XXVIII & And XXIX are self-bound dwarf 
satellites of M31, their nature and evolutionary history 
are still open questions. And XXIX seems to have 
only older (> a few Gyr) stars ()Bell et all 1201 1[ ). but 
And XXVIIFs stellar popula tion is less con strained 
by the SDSS photometry fSl ater et all 120111 ). The 
paucity of stars near And XXVIII with g — r < or 
r — i < implies it cannot have as much recent star 
formation as Leo T, but low levels of intermediate- age 
stars cannot be completely ruled out without deeper 
photometry. Further, without HI data, it cannot be 
definitively determined if these galaxies are truly passive 
dSphs, or a "transition" -type dSph/dIrr objects in a low 
star-formation episode. 

The lack of recent star formation in And XXVIII is 
surprising, as it stands in contrast to galaxies at similar 
distances from their hosts like Leo T and the Phoenix 
Dwarf. Yet unlike Tucana and Cetus, And XXVIII is un- 
ambiguously associated with M31. Thus, And XXVIII 
stands as an important data point for understanding the 
evolutionary history of dSphs satellites and their connec- 
tion to field galaxies. 
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